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Abstract The aim of the present paper is to use quasi-periodic oscillations in 
hard X-rays (HXRs) of solar flares as a diagnostic tool for investigation of im- 
pulsive electron acceleration. We have selected a number of flares which showed 
quasi-periodic oscillations in hard X-rays and their loop-top sources could be 
easily recognized in HXR images. We have considered MHD standing waves to 
explain the observed HXR oscillations. We interpret these HXR oscillations as 
being due to oscillations of magnetic traps within cusp-like magnetic structures. 
This is confirmed by a good correlation between periods of the oscillations and 
the sizes of the loop-top sources. We argue that a model of oscillating magnetic 
traps is adequate to explain the observations. During the compressions of a 
trap particles are accelerated, but during its expansions plasma, coming from 
chromospheric evaporation, fills the trap, which explains the large number of 
electrons being accelerated during a sequence of strong impulses. The advantage 
of our model of oscillating magnetic traps is that it can explain both the impulses 
of electron acceleration and quasi-periodicity of their distribution in time. 

Keywords: Flares, Energetic Particles, Impulsive Phase; Oscillations, Solar; 
X-ray Bursts, Association with Flares, Hard 

1. Introduction 

Extreme ultraviolet (EUV) images from the TRACE satellite allowed to observe 
oscillations of coronal loops during solar flares (see lAschwandenl 120041 and pa- 
pers refereed therein). Typical periods of the oscillations were 200-300 s. The 
oscillations have been explained as standing MHD oscillations of whole coronal 
loops, i.e. with nodes located at the loop footpoints. 

INakariakov et aL\ p006p observed quasi-periodic oscillations (QPOs) with 
similar periods, P = 2 — 4min, in hard X-rays (HXRs) of a solar flare and 
similar periods, P « 5 min were observed in soft and hard X-rays of a solar flare 
by IFouUon et al. \ ([2005) . In both those papers the observed periodicities were 
also interpreted as being due to the oscillations of whole coronal loops. 
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Table 1. List of limb flares witii periods P = 10— 60s^ | 



Flare 


Date; time [UT] 


P 


h 






d 


112 


B2 


No. 




[s] 


[Mm] 


[kms^-'-] 


[G] 


[Mm] 


[kms-1] 


[G] 


fb 


92/06/28; 13:54-14:00 


56; 66 


36 


4040 


184 


14.0 


785 


36 


2 


98/05/08; 01:57-01:59 


12; 18 


18 


9400 


428 


3.4 


890 


40 


3 


98/08/18; 08:18-08:23 


21; 28 


10 


3000 


136 


8.5 


1270 


57 


4 


01/05/02; 00:35-00:37 


24; (50) 


24 


6300 


290 


5.1 


670 


31 



^ Definition of parameters: P are the time intervals between HXR peaks; h is the alti- 
tude of a HXR loop-top source; v\ is the velocity estimated according to Eq. (1); B\ is 
the magnetic field strength calculated from Eq. (3); d is the diameter of a IfXR loop- 
top source; V2 is the velocity estimated according to Eq. (4); B2 is the magnetic field 
strength calculated from Eq. (5). Time interval is the interval in which the periods P 
were determined. 

^ Structure and evolution of this flare was investigated by ITomczakl II1997II . 



On the other hand, in many flares QPOs with shorter periods, P ~ 10 — 60 s, 
were observed in HXR emission (see |Lipa[ I1978|) . In the present paper we would 
hke to show how the investigation of these short-period oscillations can be used 
to investigate impulsive electron acceleration. 

In Sect. [2] observational data and their analysis are presented. Sect. [3] contains 
the discussion of our results and Sect. 0] contains our conclusions. 

In this paper we mostly investigate quasi-periodic HXR oscillations with 
periods P ~ 10-60 s, but in Sect. 12. 41 we present analysis of three limb flares 
with periods P > 120 s. 

2. Observations and their analysis 

We used HXR light-curves (recorded by Yohkoh and Compton Gamma Ray Ob- 
servatory) and HXR images (recorded by Yohkoh). First, from the Yohkoh Hard 
X-ray Telescope Catalogue ( |Sato et ai, 2006] ) we selected flares which showed at 
least three HXR impulses and the time-intervals between the impulses were of 
similar duration (see Figs. [T] and [7]). We have selected about 50 appropriate flares. 
Next, we investigated their HXR images and have taken for further analysis the 
flares for which the altitude and diameter of the HXR loop-top source could be 
reliably determined. 

2.1. Limb flares with periods P = 10 — 60s 

In this Section we investigate four limb flares. Their HXR light-curves and images 
are shown in Fig.[T] We have measured the loop-top source diameter, d, on 
various images of the same HXR source and this way we have estimated its 
measuring (r. m. s.) error to be about ±1.5]VIm. The loop-top source diameters 
and their altitudes, h, arc given in Table[TJ In the Table two values of period P 
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Figure 1. Hard X-ray observations for four limb flares. Left: HXR light-curves (dashed vertical 
lines indicate accumulation times of the HXR images shown at the right). Right: Yohkoh HXR 
images of the flares (isolines 0.12, 0.24, 0.35, 0.47, 0.59, and 0.78 Imax) and SXR images (gray 
scale). Solid line shows solar limb and dashed lines are on the solar disc. The coordinates (x, 
y) are in arcsec, the zero-point located at the disc center. 
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Figure 2. A cusp-likc (triangular) magnetic structure according to Ascliwandcn (2004). Tlie 
tiglit spiral lines show motion of trapped electrons. 




Figure 3. a Schematic diagram showing a cusp-like magnetic configuration. Magnetic fields 
PB and CP reconnect at P. The thick horizontal line is the chromosphere, b Detailed picture 
of the BPC cusp-likc magnetic structure. 



are given for each flare. They were measured as the time-intervals between the 
strong HXR peaks. The shorter period, which corresponds to higher magnetic 
field compression, is taken in the estimates in Section 2 [Eqs. (1) and (4)]. By 
multiple measuring the same time interval, P, we have estimated the r. m. s. 
errors of the measurements to be about ±2s for sharp peaks and about ±4s for 
the peaks which are not sharp or fluctuating (see error bars in Fig. [5]). 
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We have considered a hypothesis that the observed HXR oscillations are the 
result of fast kink-mode oscillations (see below about their excitation) and we 
have considered two kinds of magnetic structure where they occur: 

1. a simple magnetic loop (no inflexion points on the magnetic lines), 

2. a magnetic loop with a "cusp-like" structure at the top [Aschwanden's (|2004p 
model - see Fig. [2] . 

In the first case the observed HXR QPOs are due to oscillations of whole flaring 
loops, i.e. the nodes of standing waves are situated at the loop footpoints. 
Assuming semi-circular shape of a loop, we can calculate the velocity of the 
wave: 



The values of P, /i and vi arc given in Tablc[TJ 
For the MHD fast-mode oscillations it should be 



where va is the Alfvcn speed and Bi and p arc the mean magnetic field and 
mass density inside the oscillating loop. Hence 



IKrucker fc Linl()2008p estimated the electron number density, N, for a number 
of HXR loop-top sources. Their values were concentrated between SxlO^cm"'^ 
and 5xl0^°cm~^ with a peak around l.Ox 10^° cm~'^. Using this values of the 
electron density we have calculated values of the magnetic filed strength Bi 
which are given in Tablc[TJ 

IMariskal p006p investigated oscillations of flaring loops using Doppler-shifts 
in SXR spectra. The obtained distribution of periods showed maximum between 
3 and 5 minutes. He interpreted those oscillations as being slow-mode or fast 
magnetosonic oscillations of whole flaring loops. 

Our periods, P = 10 — 60s, arc significantly shorter than those found by 
Mariska which suggests that the oscillations seen in HXRs occur in smaller 
magnetic structures (see also comment in Section[3]). On the other hand, the 
cusp-like magnetic structure is considered to be a "standard" magnetic structure 
in solar flares (see IAschwanden| 12004^ . Therefore we have assumed that the 
oscillations seen in HXRs occur in such cusp-like structures. 

Figure[3^ shows a simple scheme of the cusp-like structure of Fig.[2l Magnetic 
field lines steeply converge around B and C. Charged particles are efficiently 
reficctcd there (magnetic mirrors). Dynamics of the magnetic field inside the 
cusp volume BPC is explained in Fig.jSb. Magnetic field PB and CP reconnect 
at P. This generates a sequence of magnetic traps which move downwards and 
collides with a stronger magnetic field R below. The collision causes compression, 
the magnetic and gas pressure increases inside the traps, so that the compression 
is stopped, the traps can expand and perform magnetosonic oscillations. 

The excited fast magnetosonic waves propagate along the magnetic field lines, 
they meet regions of steep increase of Alfven speed near B and C, where they 




(2) 




(3) 
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are efficiently reflected back into the trap. Hence, the magnetosonic waves will 
be efficiently trapped within the magnetic traps and we can describe them as 
standing fast magnetosonic (kink) oscillations trapped between B and C. 
We can calculate the velocity of the waves 

V2 = 2s/P, (4) 

where s is the distance between B and C taken along the magnetic field lines. 
Diameters, cZ, of the HXR loop-top sources have been used to estimate values of 
s: 2s ~ ird. Hence 

nd , , 

V2 = -^ (5) 

Next we have calculated expected values of the magnetic field, 

(6) 

The values of V2 and B2 are given in Table[T] The values show surprisingly low 
dispersion (see also values in Table[2]). We consider this to be a confirmation that 
the oscillation are concentrated in cusp-like structures. 

This model well describes generation of the short-duration HXR impulses, see 
IJakimiecl (|2002i 12005^ . During the collision of the downwards moving magnetic 
field with the magnetic field below, the magnetic traps seen in Fig. [5)3 undergo 
compression and the magnetic mirrors (Mi , M2 in the Figure) move towards the 
center of the trap. Particles contained in the traps are accelerated by betatron 
effect and Fermi acceleration due to collisions with the converging magnetic 
mirrors. Let us denote the trap ratio x = Bmax/ Bmim where B„iax is the 
magnetic field strength at the mirrors B and C and Bmin is the strength at 
the middle of a trap. At the beginning of compression values of x are high, they 
decrease during the compression and accelerated particles can escape towards the 
loop footpoints where they generate the HXR footpoint emission. Some electrons 
emit HXRs before they escape from the trap and this is why we can see the 
cusp-like volume BPC as a HXR loop-top source. The parameter x reaches its 
minimum value, Xmim during the maximum of compression. The values Xmin 
can be different for different traps seen in Fig. [3)3 and the traps with low values of 
Xrnin wiU bc most efficient in particle acceleration and precipitation. Traps which 
have higher values of Xmin are most efficient in maintaining the oscillations. 

Numerical simulations of this model of particle acceleration were carried out 
by |Karlicky fc Kosugi (|2004p . A similar model of particle acceleration was 



proposed by Somov & Kosugi (|1997p . The main difference is that they assumed 
that the downward reconnection flow is super- Alfvenic and therefore a fast shock 
is formed near the bottom of the cusp-like volume. However, observations of 
lAsai et all (|2004p show that the reconnection downflows can be sub- Alfvenic. 
Our model of particle acceleration is valid also in the case of sub- Alfvenic ve- 



locities in the cusp-like volume. [Bogachev &: Somov] (|2005p have carried out 



analytical calculations of particle acceleration in collapsing magnetic traps. 



SOLA: jj+mt.sp.tex; 5 August 2009; 12:40; p. 6 



QPO in HXRs of solar flares 



Table 2. List of flares located on the solar disc" ■ 



Flare 


Date; time [UT] 


P 


D 




Bi 


d 


1>2 


B2 


No. 




[s] 


[Mm] 


[kms^-'-] 


[G] 


[Mm] 


[km s^-"-] 


[G] 


5 


91/10/31; 09:09-09:11 


30; 39 


21 


2160 


99 


5.5 


574 


26 


6^ 


93/05/14; 22:04-22:05 


25; 33 


29 


3600 


166 


5.3 


660 


30 


7 


93/06/07; 14:20-14:22 


60; 68 


32 


1650 


76 


13.7 


717 


33 


8 


94/01/16; 23:16-23:18 


51; 51 


19 


1180 


54 


12.3 


757 


35 


9 


99/07/28; 08:13-08:14 


31; 35 


16 


1660 


76 


6.8 


690 


32 


10 


00/06/23; 14:25-14:26 


35; (76) 


17 


1530 


70 


6.8 


610 


28 


11 


00/11/25; 18:38-18:39 


21; 28 


22 


3270 


150 


4.1 


613 


28 



" Definition of parameters: D is the distance between the HXR footpoints; vi = ivD/P. 
The other parameters are the same as in Tablc[T] 

^ Structure and evolution of this flare was investigated by IKolomahskil l)2007|l . 




Figure 4. Observations of the flare of 16 January 1994. a GRO/BATSE light curve (nominally 
25-50 keV); b YoMoh HXT image, 23-33 keV (isolines) and SXT image (gray scale). 
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2.2. Flares located on the solar disk 

Usually, HXR footpoints sources are well seen in the case of flares located on the 
solar disk, ft is more difficult to recognize HXR loop-top sources and to estimate 
their altitudes. 

We have selected a sample of disk flares for which QPOs were observed and 
loop-top sources could be identified and measured (see example in Fig. 3]). The 
distance, D, between the HXR footpoints has been taken for these flares as the 
measure of the loop size {h = D/2 for a semicircular loop). We have calculated 
wave velocity, wi, for a whole loop as 

vi = ttD/P (7) 

and the corresponding magnetic field, Bi, from Eq. (3). 

Next, we have calculated the values of V2 and B2 for the cusp-like model from 
Eqs. (4) and (5). The obtained values are given in Tabled 

We see that Tabled] confirms the properties seen in Table[TJ The values of vi 
and Bi show large dispersion. The values of V2 and B2 show much smaller dis- 
persion and they are in very good agreement with those in Tablc[T] We consider 
this to indicate the advantage of the cusp-like model. 

2.3. Investigation of sequences of HXR images 

Important information can be derived from analysis of sequences of HXR images. 
Such a sequence for a flare of 16 January 1994 is shown in Fig. [5] {Yohkoh 23- 
33keV images). The location and structure of this flare are shown in Fig.|4l The 
loop-top (LT) and the footpoint (Fl and F2) sources are well seen in the images 
of Fig.O At the top of this figure is shown the HXR light curve with the time 
intervals for which the images Nos. 1-6 have been accumulated. Images Nos. 2, 4, 
and 6 correspond to the main HXR peaks and images Nos. 3 and 5 correspond 
to "valleys" between the peaks. 

Let us note that - in spite of moderate spatial resolution - a trace of a 
triangular structure of the HXR loop-top source is seen in images Nos. 1 and 3 
of the figure. It is not seen during the maximal compression of the traps (images 
Nos. 4 and 6). 

Figure[B] shows time- variation of the intensity of the brightest pixel in foot- 
points Fl and F2 and in loop-top source. We see that the main HXR peaks, 
Nos. 2, 4 and 6, are best seen in the footpoint intensity, but they are less marked 
in the loop-top intensity. This is in agreement with the outlined model: high 
energies of electrons are achieved at the end of trap compression and then 
they can easily escape towards the footpoints, so that their contribution to the 
loop-top emission is moderate. 

Let us note that there is some asymmetry in time-variation of the footpoint 
intensity (Ji > I2 at peak No. 4 and Ii < I2 at peak No. 6). This indicates that 
there is some asymmetry in the oscillation of main trap (left-right oscillations in 
Fig. [3)3) which causes that the time variation of the ratio x is somewhat different 
for the two ends of the trap. 
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Figure 6. Time-variation of the HXR intensity (23-33 keV) of the brightest pixel in the 
footpoints Fl and F2 and in the loop-top source. 



Similar results for other flares will be presented in another paper. 

Hitherto we have considered impulses (peaks) of HXR emission. But it is 
important that HXR emission from the loop-top and footpoint sources is seen 
also during the HXR valleys (see Figs.[5][6l). This means that the process of elec- 
tron acceleration and precipitation is operating all the time during the impulsive 
phase, also during the HXR valleys. To explain this, let us return to Fig.[5]D. Be- 
ginning of intensive reconnection at P generates vigorous reconnection downflow 
which collides with the transverse magnetic flcld R and causes compression of 
a "main" magnetic trap and onset of its oscillations. But the reconnection is 
continued and this means that a sequence of magnetic traps, which are located 
above the main trap, undergo compression and can oscillate. The compression 
and oscillations of these traps are shifted in time relative to main trap, so that the 
acceleration and precipitation of the electrons occurs during the whole impulsive 
phase, including the HXR valleys. Hence, the HXR observations suggest that the 
cusp-like volume BPC is a dynamic region which is filled with moving, oscillating 
and colliding magnetic traps and it is the volume where electrons are accelerated. 

If the main oscillating trap is a dominant one (i.e. it contains large number of 
particles), quasi-periodic sequence of HXR impulses will be observed. If, however, 
there are several oscillating traps of similar contribution to the HXR emission, 
the HXR peaks can be randomly (i.e. not quasi-periodically) distributed, as it 
is the case for impulsive phase of many flares. 
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ITomczakl (|200ip investigated positions of HXR and SXR loop-top sources 
for limb flares, during their impulsive phase. He has found that in most cases 
the HXR and SXR sources are co-spatial or nearly cospatial. This result has 
been recently confirmed by IKrucker fc LinI (|2008p using RHESSI observations; 
see also ITomczakl p009p . These results indicate that in most cases both the 
HXR and SXR loop-top emission, recorded during the impulsive phase, comes 
from the same or nearly the same plasma volume. Let us note that this common 
volume of loop-top HXR and SXR sources confirms that HXR loop-top sources 
(like SXR ones) are also feeded by chromospheric evaporation upflows. 

iMasuda, fil994|) has found that in some flares the HXR loop-top source is 
located clearly above the SXR loop - see his flare of 13 January 1992 at 17:28 
UT. The schematic diagram in Fig. [2] just concerns the "Masuda flare" of 13 
January 1992. 

In order to explain such cases as the Masuda flare we should take into account 
that in the BPC cusp volume there is a sequence of many moving traps and 
the dynamics of compression of individual traps is determined by the time 
variation of the velocity, v{t), of reconnection outflow at the top of the cusp. 
The traps in the lower part of the cusp volume are "typical" : they undergo 
efficient compression, i.e. Xmin is low for these traps. HXR emission from these 
low traps was seen in the Yohkoh L-channel (14-23 keV) - see IMasudal (fT994l) 
. But these traps were not seen at higher energies {E > 20keV) because (1) 
the higher energies of electrons are achieved at the end of trap compression and 
then they can immediately precipitate toward the loop footpoints, (2) there was 
a dominant HXR {E > 20keV) sources at higher altitude. The heating of the 
footpoints causes filling of the loop with hot plasma and hence the loop is seen 
in SXRs. This description is confirmed by the fact that HXR footpoint sources 



were located just at the ends of the SXR loop (see Fig. 12.12 in IAschwanden| 




The enhanced HXR emission {E > 20keV) from the upper part of the cusp 
indicates that the compression of the traps located there is less efficient, so that 
the ratio x remains high also during maximum of compression and hence the 
escape of accelerated electrons is more difficult (the traps are "semiclosed" ). The 
accelerated electrons remain within the trap for longer time and they emit HXRs 
before they escape from the trap. 

Hence, the main point of our interpretation of the Masuda flare is that the 
footpoint HXR sources are fed by the accelerated electrons coming from "typical" 
traps, which are situated low in the cusp volume BPC, and the top HXR source 
is a "semiclosed" magnetic trap, from which escape of particles is more difficult. 

Our model can also easily explain the cases when HXR loop-top source is 
located below the SXR loop. Then "semiclosed" traps are located at the bottom 
of the BPC cusp volume, below the "typical" traps which are seen in SXRs. Such 
examples have been found by IKrucker fc LinI (|2008p and ITomczakl (|2009p . In 
their papers Ah is the vertical distance between the centroids of HXR and SXR 
loop-top sources and Ah < means that the HXR source is located below the 
SXR one (see Fig. 2 in IKrucker fc Linll2008p . Let us note that such cases cannot 
be easily explained by models which assume that a fast shock is formed within 
the BPC cusp volume l|Masuda| [T9941 [S omov fc Kosugi[ I1997p . 




mm- 
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Figure 7. Hard X-ray observations for three limb flares. Left: HXR light-curves (dashed 
vertical lines indicate accumulation times of the HXR images shown on the right). Right: 
Yohkoh HXR images of the flares. Solid line shows solar limb and dashed lines are on the solar 
disc. 
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Table 3. Limb flares with periods P > 120 s'' | 



Flare 


Date; time [UT] 


P 


h 


"1 


Bi 


a 


V2 


B2 


No. 




[s] 


[Mm] 


[kms^^] 


[G] 


[Mm] 


[kms"-"-] 


[G] 


12 


93/03/02; 21:10-21:15 


(93); 137 


38 


1700 


78 


15 


690 


32 


IS'' 


98/09/30; 13:29-13:38 


140; 166 


38 


1700 


78 


20 


900 


41 


14b 


01/04/03; 08:18-08:23 


150; 206 


34 


1400 


64 


20 


840 


39 



Definition of parameters: a is the size of cusp-hkc structure; V2 = 2na/P; the other 
parameters are the same as in Table[T] 

^ Structure and evolution of these flares were investigated by IBak-Steslickal I I2007II . 



2.4. Limb flares with periods P > 120 s 

In this Section we consider three hmb flares with periods P > 120 s (Fig. [7] and 
Table[3|). AU they are high structures (see the altitudes, h, in Table[3]) and long- 
duration events (LDE). Yohkoh SXR observations were available for the flares 
Nos. 13 and 14. Their SXR images (grey scale) together with HXR isophotes are 
shown in Fig. [SI The SXR images show that they were big arcade flares. 

In particular, very interesting is the large triangular structure seen in the HXR 
image of the flare No. 13 of 30 September 1998 (see Fig.[7|). The dominant HXR 
sources, B and C, are located at the altitude of SXR arcade (see Fig. [8]). Weak 
HXR footpoint sources, located on the solar disc, are better seen in Yohkoh Ml 
(23-33 keV) and M2 (33-53 keV) images. 

We have assumed that this triangular structure is analogous to BPC cusp- 
like structure in Fig. [3] where the observed oscillations occur. The two dominant 
HXR sources, B and C, correspond to the magnetic mirrors, B and C, in Fig.[3l 
Charged particles spend somewhat longer time (calculated per unit length along 
the magnetic trap) near the magnetic mirrors than in other parts of the trap and 
therefore the HXR emission is enhanced near the mirrors. The enhanced emission 
near P and along the line running downwards from P, is due to the fact that 
plasma density is somewhat higher there because of reconnection stream flowing 
from the current sheet at P. The higher density increases the HXR emissivity. 

We have taken the distance, a, between the top source P and the line con- 
necting B and C, as the height of the cusp structure, and calculated the velocity 
V2 as 

= ^ (8) 

and then the magnetic field strength B2 from Eq. (5). The values of vi and Bi 
calculated from Eqs. (1) and (3), as well as V2 and B2 are given in Table[3l 

Flare No. 14 of 3 April 2001 is a similar case. It was also a large arcade flare 
and the sources B and C were located at the altitude of SXR arcade. A top 
HXR source P is seen in Yohkoh Ml (23-33 keV) images which conflrms that 
this is also a similar triangular HXR structure. HXR footpoint sources are seen 
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Figure 8. Soft X-ray images of the flares Nos. 13 and 14 (gray scale). The isophotes show 
HXR images Hke in Fig. [7] it is 14-23 keV image in the case of 30 September 1998 flare and 
23-33 keV image for the 3 April 2001 flare. 
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on the solar disc, but they are weaker than the B and C sources. Again, we have 
measured the height of the cusp structure as the distance between the top source 
P and the hne connecting B and C and the obtained values of vi, V2, and 
B2 are given in Table[3l 

Our model (connection of sources B and C by a cusp-like structure like in 
Fig-Eb) explains why in flares Nos. 13 and 14 the sources B and C are similar 
in intensity and their time evolution is similar. Figure[8] indicates that in the 
both flares the plane BPC of the cusp-like structure is not perpendicular to the 
arcade axis. This can be due to the shear in footpoint motion at the photosphere 
(see Fig. 3.10 in IPriest"! I1987p . The sources B and C arc located at the places 
where cusp-like magnetic lines meet the arcade loops (the places are located on 
the opposite sides of the arcade axis). 

No SXR images were available for the flare No. 12 of 2 March 1993. In the 
HXR image of this flare (Fig.[7|) strong footpoint sources are seen on the solar 
disc and a triangular structure BPC is seen at the top. The strong footpoint 
emission indicates that the oscillations of magnetic traps were more dynamic 
(low value of Xmm) and therefore the accelerated electrons could easily escape 
towards the footpoints. We have measured the altitude, h, and the height, a, of 
the cusp-like structure at the top and obtained values of vi, Bi, V2, and B2 are 
given in Table[3l 

In Table[3] we see that the values of f 2 and B2 are similar to those obtained 
for flares with periods 10-60 s (Tables[l] and [2|) . This supports our assumption 
that the large triangular HXR structures, considered in this Section, are large 
magnetic cusp-like structures. 

3. Discussion 

A lower limit for the Alfven speed can be estimated as follows: 

= 2.2 X 10"S/ViV [cms-i] (9) 

where N is the electron density. For a magnetic trap it should be 

B^/{8tt) > 2NkT, (10) 

where k is the Boltzmann constant and T is the temperature. This gives 

B/Vn> 8.3 X 10^ Vt (11) 

Hence the lower limit for the Alfven speed is [Eqs. (9) and (11)]: 

VA > 1.8 X W^Vt [cms-i] (12) 

Taking T > 10^ K for the temperature in the trap, we obtain 

WA> 5.7x10^ cms-i. (13) 

We see in Tablcs[TH3]that our estimates of the Alfven speed, V2, are higher than 
this lower limit, but they are not much higher. This means that the magnetic 



SOLA: j j+mt_sp.tex; 5 August 2009; 12:40; p. 15 



J. Jakimiec & M. Tomczak 



50 



40 r 



30 r 



20 r 



10 r 



L 



+ 

+ + 
+ + + 



50 100 150 

P [sec] 



Figure 9. Correlation between the HXR loop-top source diameter, d, and the oscillation 
period, P. For large sources [P > 120 s) the values of 2a were taken for d - see text. Point 
marked with a square relates to the flare investigated by Ofman & Sui (2006). The larger 
cross at the first data point on the left shows r. m. s. error bars which concern all small HXR 
sources, P < 60s, marked with small crosses. 



pressure, B^/(87r), is not much higher than the gas pressure, p, in the acceleration 
volume. This is important for MHD modeling of flaring loops. 

In Fig. [9] the diameter, d, of the HXR loop-top source is plotted vs. the os- 
cillation period, P. For short periods, P = 10 — 60 s, good correlation between 
d and P is seen. [A diagram displaying h (or D/2) vs. P shows much larger 
dispersion.] This is a confirmation that the investigated oscillations occur in 
confined magnetic structures displayed as the HXR sources. This supports our 
model of the oscillations of magnetic traps in cusp-like magnetic structures. 

For the large HXR sources (P > 120 s), the values of 2a were taken instead of 
d in Fig.[9l where a is the height of a triangular HXR structure. It appears that 
the large HXR structures also fit the correlation between d and P determined 
by the small structures. This supports our assumption that the large triangular 
HXR structures are the cusp-like magnetic structures in which the oscillations 
with period P > 120 s occur. 

lOfman fc Suil (|2006p have analyzed RHESSI observations of QPOs in a flare 
of 19 January 2005 with period P ~ 150 s. We have measured the size, a, of 
the HXR loop-top source in the same manner as for our three large Yohkoh 
flares and we have plotted additional point marked with a square in our test 
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diagram in Fig.[9l We see that this additional flare fits quite well into our model 
of oscillating magnetic traps. 

It may be important that all four large flares in Fig.[5]have periods about 140- 
150 s, which suggests that they may be overtones of the 300 s solar oscillations. 
This would mean that the oscillations are maintained due to a resonance with 
the 300 s oscillations. In particular, this seems to be probable in the case of 
our flares Nos. 13 and 14 and the flare investigated by lOfman fc Suil (|2006p 
, in which the HXR oscillations are quasi-sinusoidal (not typical sequences of 
increasing impulses like in Fig. llOp . 

For a number of flares IMariskal l|2006p has found decaying oscillations of 
periods 5.5±2.7 minutes in Doppler shifts of SXR spectral lines. He explained 
them as being due to slow-mode or fast magnetosonic oscillations of whole flar- 
ing loops. We have investigated Yohkoh HXR records for those flares and we 
have found no counterpart of the Mariska's oscillations in the HXRs. Mariska's 
flares were rather small, hence the observed oscillations most probably were 
the relaxation oscillations of whole flare magnetic structure which were excited 
by collision of violent chromospheric evaporation upflows. Moreover, for some 
Mariska's flares we have found HXR quasi-periodic oscillations of period 10-60 s 
analogous to those investigated in the present paper. This is an additional indi- 
cation that the HXR oscillations with periods 10-60 s occur in smaller volumes 
and they are different from the longer-period oscillations observed by Mariska. 
These cases will be investigated in a separate paper. 

Important information can be obtained from CGRO/BATSE observations 
due to their high sensitivity in HXRs. An example of such observations is shown 
in Fig.U] (upper panel). We see that a sequence of weak quasi-periodic impulses 
began about 23:09 UT, but a typical impulsive phase began about 23:14UT. 
(Similar weak oscillations before the impulsive phase are seen in flare of 2 March 
1993 in Fig. [7]). These observations are in agreement with the model of oscillat- 
ing traps: The oscillations began about 23:09 UT. Accelerated particles hit the 
footpoints and chromospheric evaporation begin. 

The plasma moves along the same magnetic lines, connecting traps with foot- 
points, and about 23:14UT it reaches the oscillating traps. During the expansion 
of a trap this plasma fills the trap and therefore a larger number of electrons 
will be accelerated during the next compression of the trap. Hence, there is now 
a feedback between the number of accelerated electrons and the density of the 
plasma which reaches the traps. This feedback causes that the HXR intensity 
increases with time, as it is seen in Fig. 2] 

Hence, the model of oscillating magnetic traps provides a simple explanation 
of the well-known paradox that the number of electrons, which arc accelerated 
during big flares, is much higher than it can be provided by surrounding corona. 

Let us repeat that, according to discussion of Fig. 21 chromospheric evapo- 
ration begins a few minutes before the beginning of the impulsive phase. The 
beginning of the evaporation causes increase of the SXR emission. Hence, the 
model of electron acceleration in oscillating magnetic traps is in agreement with 
the well-known fact that SXR emission of a flare usually begins a few minutes 
before the beginning of HXR impulsive phase. 



SOLA: j j+mt_sp.tex; 5 August 2009; 12:40; p. 17 



J. Jakimiec & M. Tomczak 



GRO/BATSE DISCLA Rates (1.024 sec averages) 
3.5x1 o"^ I _ 




21:58 22:00 22:02 22:04 22:06 22:08 



Figure 10. HXR light curve (nominally 25-50 kcV) for the flare of 14 May 1993 (flare No. 6). 

Thermal conduction also contributes to the energy transport from the loop 
top to its footpoints during the impulsive phase and it increases the rate of 
chromospheric evaporation. 

Our model of oscillating magnetic traps explains very well why the time 
intervals, P, between the HXR impulses are not equal. The oscillations occur 
in a volume of changing magnetic field and plasma density, so that the Alfven 
speed can change from one pulse to another. 

In the present paper we have considered a cusp-like magnetic configuration. 
We see from this example that two factors are necessary for magnetic-trap 
compression and excitation of their oscillations: a magnetic-field reconnection 
and an obstacle (a transverse magnetic field) on the way of the reconnection 
flow. But the location and orientation of the reconnection site can be different 
in different flares. 

Let us note that our model of exciting oscillations of magnetic traps is similar 
to that proposed for oscillations during onset of magnetospheric substorms. The 
substorm begins with a violent reconnection in the magnetospheric tail (see 
observations presented by INakamural I2004p . The reconnection outflow directed 
Earthward collides with the inner (dipole-like) magnetic field, causes magnetic 
field compression and excitation of oscillations. An example of magnetic field 
oscillations during the onset of a substorm is shown in Fig.[TT](a smooth substorm 
variation of the magnetic field has been removed from these records). We see 
a short sequence of oscillations which is similar to what we see in HXRs of 
investigated flares (the smooth variation of HXR intensity was not removed 
from our light curves). 
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Figure 11. Magnetic field oscillations during the onset of a geomagnetic substorm of 23 
November 1970 (smooth variation has been removed) as recorded by seven geomagnetic stations 
in Canada. The distance between the horizontal lines amounts to 316 7. From Nishida (1978). 

4. Summary and conclusions 

We have investigated flares for which quasi-periodic HXR oscillations and HXR 
loop-top sources were clearly seen. Our analysis has shown that the observed 
HXR oscillations are confined within the loop-top sources. This is supported by 
a correlation between the periods, P, and the sizes, d, of the loop-top sources. 

IMariskal (^005)) has found that oscillat ions of whole flaring loops had periods 
of 5.5±2.7 minutes. Recently we have found that HXR oscillations with periods 
P ~ 10 — 60 s occurred in some of the Mariska's flares which confirms that the 
oscillations seen in HXRs, occurred in smaller magnetic structures (these results 
will be presented in a separate paper). 

We have argued that a model of oscillating magnetic traps is adequate to 
describe the HXR observations. During the compression of a trap particles 
are accelerated, but during its expansion plasma coming from chromospheric 
evaporation, fills the trap. This feedback between the particle acceleration and 
filling of the traps explains high number of electrons which are accelerated in 
strong flares. This model is clearly supported by finding weak HXR oscillations 
a few minutes before the impulsive phase (see Sect.[3|). 
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We have estimated the mean strength of the magnetic field inside the oscil- 
lating traps, ft amounts to about 30-40 Gauss, which is lower than assumed in 
other papers (see lAschwandeni I2004p . 

We have compared the properties of the ffXR oscillations with available mod- 
els of flare magnetic structure and particle acceleration and indicated the most 
adequate scenario (see Sections [3] and • 

fIXR light curves can be divided into two components (see lAschwanden et al.[ 
If 996| ITomczaicl I200ip : a smooth component (a smooth line connecting valleys) 
and impulsive one (increases above the smooth line). The smooth component is 
the result of continuous acceleration in many traps whose oscillations are shifted 
in phase (this is due to continuous reconnection at the top of cusp structure). The 
increase of the smooth component during the impulsive phase (like in Fig. fTO]) 
is parallel to increase of the SXR loop-top emission which confirms that the 
increase of the smooth HXR component is also due to the increase of the density 
within the cusp-like volume. 

If there is a pronounced maximum in the time-profile of the reconnection rate 
at the top of cusp, then a dominant (large) trap will be generated within the 
cusp-like volume and its oscillations will be seen as quasi-periodic oscillations in 
HXRs. 

We have pointed out that in most cases the HXR and SXR sources are cospa- 
tial or overlap during the impulsive phase which confirms that both emissions 
come from the common cusp-like plasma volume. But in seldom cases the HXR 
and SXR loop-top sources are clearly separated (like in "Masuda flare"). We 
explain such cases of HXR sources as being due to magnetic traps, for which 
the ratio x remains high during the compression, so that accelerated electrons 
remain there for a longer time and therefore generate enhanced HXR emission. 
This explanation is also important for the flares having weak HXR footpoint 



emission (like in the case of Veronig &: Brown} I2004p 



In our paper of 1998 ( [Jakimiec et ai, 1998[ ) we have found that (a) loop- 
top flare kernels are multithermal, (b) significant random motions occur in the 
kernels. On this basis we have proposed that the kernels are filled with MHD 
turbulence. In the present paper we have used observations of HXR oscillations 
to improve the model of loop-top kernels. We have come to the conclusion that 
during the impulsive phase the kernels are filled with oscillating magnetic traps. 
The oscillations are excited and maintained by violent reconnection flow. The 
advantages of this model arc the following: (f ) it incorporates a simple model 
of particle acceleration (Fermi and betatron acceleration), (2) the particles can 
easily escape after their acceleration. We will observe this plasma volume as being 
"turbulent" (random motions due to oscillations shifted in phase; temperature 
differences between different traps) . (We think that MHD turbulence can develop 
later on due to collision of chromospheric evaporation upflows. This would ensure 
long duration of the kernels.) 
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